Critical phenomena in disc-percolation model and its application to 

relativist ic heavy ion collisions 
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Through studying the critical phenomena in continuum-percolation of discs, we find a new ap- 
proach to locate the critical point, i.e. using the inflection point of Poo as an evaluation of the 
percolation threshold. The susceptibility, defined as the derivative of Poo, possess finite-size scaling 
property, where the scaling exponent is the reciprocal of v — the critical exponent of correlation 
length. The possible application of this approach to the study of the critical phenomena in rela- 
tivistic heavy ion collisions is discussed. The critical point for deconfinement can be extracted by 
the inflection point of -Pqgp — the probability for the event with QGP formation. The finite-size 
scaling of its derivative can give the critical exponent v, which is a rare case that can provide an 
experimental measure of a critical exponent in heavy ion collisions. 
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I. INTRODUCTION 

The general feature of the phase diagram of strongly 
interacting matter has become increasingly well estab- 
lished Following a region of crossover around a tem- 
perature of 170 - 190 MeV at zero baryon chemical po- 
tential fig, increasing [i^ leads to a critical point, beyond 
which the system shows a first order transition from con- 
fined to deconfined phase. 

Recently, the interest in QCD critical point progres- 
sively rises. With the plan of RHIC low energy scan and 
GSI new facility, the study of relativistic heavy ion colli- 
sions becomes concentrating on searching for the critical 
point and observing the relevant critical phenomena, in 
particular measuring the critical exponents. 

On the theoretical side, studying the critical phenom- 
ena in chiral symmetry restoration is available [3L Many 
discussions about universality have been active [3|. The 
chiral condensate and Polyakov loop are proposed as or- 
der parameters for chiral restoration and color decon- 
finement, respectively. However, it is a pity that both 
of these variables can not be directly measured by ex- 
periment. What observable should we measure, how to 
locate the critical point and extract the corresponding 
critical exponents in heavy ion collisions have not been 
clear so far. 

On the experimental side, what we observed in heavy 
ion collisions is the deconfined partonic degree of free- 
dom Q. Despite the fact that confinement is a long 
standing problem not solved by theory, studying the criti- 
cal phenomena in deconfinement phase transition is more 
realistic for experiments. Since chiral symmetry is hard 
to measure and deconfinement is easier to observe, we 
propose to phenomenologically study the critical phe- 
nomena from the deconfinement aspect and give some 
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hints for experiments. 

Principally speaking, it is impossible to get the crit- 
ical point because of the limited system size in rela- 
tivistic heavy ion collisions. In this paper, by means of 
studying the critical phenomena in finite-size continuum- 
percolation of disks, we find that P^, the probability for 
event with infinite cluster, has an inflection point, which 
is a good approximation for the critical point. Finite- 
size scaling method is further used to extract the crit- 
ical exponents v from the distribution of susceptibility. 
Since percolation has some resemblance with deconfine- 
ment [5j, a possible application of the method used in 
percolation to the study of the critical phenomena in rel- 
ativistic heavy ion collisions is discussed. 



II. THE CONTINUUM-PERCOLATION OF 
DISC 

The continuum-percolation problem has various for- 
mulations, among which the problem of spheres is most 
popular. Equally sized spheres are placed at random 
in a substrate. The spheres support the transport and 
the substrate does not. When we put enough spheres 
in the substrate, the overlapping spheres form an infi- 
nite cluster, and the system is able to support a long- 
range current. This model has been well studied and was 
used to describe hopping conduction in doped semicon- 
ductors @ and phase transitions in ferromagnetics 0]. 
Such a model together with its critical behavior also bare 
some resemblance with the transition from hadron gas to 
quark-gluon plasma. 

In QCD, hadrons are color-singlet bound states of more 
basic colored objects — quarks and gluons. Hadronic 
matter, consisting of colorless constituents of hadronic 
dimensions, can turn at high temperature and/or den- 
sity to a quark-gluon plasma of colored quarks and gluons 
as constituents in a much larger volume. This deconfine- 
ment transition leads to a color-conducting state and thus 
is the QCD counterpart of the insulator-conductor transi- 
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tion in atomic matter. Suppose hadrons have an intrinsic 
size. At low density, we have a hadron gas. When this 
becomes so dense that the formation of an infinite clus- 
ter occurs, it turns to a quark-gluon plasma. In this case, 
the connectivity (cluster formation) determines the dif- 
ferent states of many-body systems. The lesson learned 
from spin systems also indicates that cluster formation 
and associated critical behavior are the more general fea- 
ture [f|. 



A. Model parameters 
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FIG. 1: (Color online) for different system size. From left 
to right the system sizes are 25, 100, 300, respectively. The 
error bar shown is systematic error induced by the variation 
of parameter 9q. Statistical errors are small and not shown. 
The dashed curves are the susceptibility defined in Eq. (4) 



In relativistic heavy ion collisions the incoming nuclei 
become discs of vanishing thickness due to Lorentz con- 
traction and the collision region is two-dimensional. So 
we turn to study two-dimensional percolation. 

In a two-dimensional system, spheres become discs. 
Discs with radius a, called cells, distribute randomly in 
the system of a big disc with radius R. The ratio R/a 
determines the system size and is denoted by s. The con- 
trol variable for percolation is ry — the ratio of the total 
area of all the cells and the area of the big disc, i.e. 

Ntto 2 No 2 

where N is the number of cells. The increase of cell 
number N leads to the increase of r\. The value of rj at 
which an infinite cluster appears is the critical point and 
is denoted by 7y c . 

The definition for infinite cluster varies. In finite sys- 
tem a cluster spanning the system is called an infinite 
cluster. Let (r,,^,) be the polar coordinates of cell i. If 
ri <E [R — a, R] we call the cell i as a boundary cell. If 
cells i and j are two boundary cells belonging to the same 
cluster, and \9i — 6j\ > 6q then we say the current cluster 
is infinite. Here 8q is a parameter. 

Thus the control variable is r\ and the model parame- 
ters are s and 9q. 



B. Locating the critical point 

When s and 8q take certain values respectively, the 
probability for the event with infinite cluster, denoted by 
Poo, is an increasing function of r). For infinite system 
size, Poo = holds for all r\ <r\ c . When r\ > r) c , Poo picks 
up and sharply increases to unity. System of infinite size 
is unrealistic. Poo for finite system of three different sizes 
are measured and shown in Fig. [TJ 

Different parameters, e.g. 8q = 90°, 135° and 175°, 
are tried and the induced systematic errors are shown as 
error-bars in the figures. The solid curves in Fig. [1] are 
the fitting result by using the function 

Poo( ? y)= 1 + taDh[ : i(?? ' C2)] . (2) 
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FIG. 2: (Color online) The behavior of C2 versus system size. 
A tendency to saturation can be seen. 



It can be seen that for all three cases the fitting function 
always lies within the systematic errors. The fitting func- 
tion, i.e. Eq.([2l), has an inflection point at 77 = C2. As s 
increases, P^ tends to a step function and the inflection 
point C2 tends to r\ c . For finite system size, the inflection 
point "q = C2 can be used as an evaluation of i] c . 

The behavior of C2 versus different system size is shown 
in Fig. [5] A tendency to saturation can be seen. The sat- 
uration value of C2 is 1.1198 ±0.0047, i.e. the percolation 
threshold is r] c = 1.1198 ± 0.0047. 

For a nucleus-nucleus collision, the contraction induced 
by relativistic motion makes the colliding nucleus like a 
disc, in which disc-like nucleons percolate. For central 
Au-Au collision, the radius of the interacting region is 
about 7 fm and the hard core of each nucleoli is about 
0.1 fm, resulting in a system size of about 70. The in- 
flection point for s = 70 is about 1.1054. Using it as 
an approximation for r] c , the error is 1.3%. So using the 
inflection point as an approximation for critical point is 
easy to measure and exact enough. 

In heavy ion collisions, when hadrons connect to form 
infinite cluster, color conducts in the system and quark 
gluon plasma is formed. If the formation of QGP can 
be signaled in each collision, then the probability for the 
event with infinite cluster, denoted by Poo, becomes the 
probability for the collision with QGP — — Pqgp • Pqgp 
should be measurable in experiment, provided that an 
available signal is found, and thus in principle the critical 
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FIG. 3: (Color online) P^ sual (open circles) as a function of 
n with system size s = 1000. The solid curve is the fit of the 
lower part to Eq. ([3]). For comparison the P x for s — 1000 is 
also shown as triangles . 



FIG. 4: (Color online) The distribution of susceptibility n at 
critical point r\ c as a function of system size s. The solid line 
is a fitting by Eq. ©. 



point can be obtained according to the inflection point 
ofP QG p. 

Poo is also referred to by some authors [8[ as percola- 
tion cumulant, in analogous to the Binder cumulant 
which will intersect at the critical point for different sys- 
tem sizes. Extracting the intersection point of Poo is a 
method to estimate the critical point. However, to use 
the inflection-point method proposed above only needs 
the collisions of one kind of ion at various energies, while 
to use the intcrsecting-point method need the collision at 
different energies of more than one kind of ion. There- 
fore, the inflection-point method is more realistic and 
the value of r\ c obtained for system sizes as small as Pb 
(s = 71), Au (s = 70) is a good estimation with error 
about 1.3%. 

The inflection-point method is also applicable for the 
P-infinity defined in another way. In Ref. [l0( P-infinity is 
defined as the probability of a cell belonging to an infinite 
cluster, which will be denoted by P^ ua[ in the following. 
Fig. [3] shows P^ sual as a function of rj when s = 1000. In 
case of s = 1000, the system is a big one and Poo rises like 
a step function, cf. Fig. [3l but P^ ued has a step function 
sharp rise only at the bottom side (Po^ 8ual = 0), while at 
the upper side Po^ sual tends to unity smoothly. 

In this case, we could only fit the left part of P t i c 3ual ( 7 ?) 
to Eq. ([2]), or alternatively, to a modified fitting function 
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ci [c 2 + tanh (c 3 {q - c 4 ))] 



(3) 



The fitting is shown as solid curve in Fig. [31 The inflec- 
tion point of this half-fitted curve can also be used as an 
approximation of the critical point. 



C. Susceptibility and the critical exponent v 

Susceptibility is defined by the response of the system 
to small external forces. We define the susceptibility in 
percolation as the derivative of Poo with respect to the 
control variable 77, i.e. 



K{r),s) 



dPoo{ri,s) 
dr/ 



(4) 



k shows a peak at the inflection point of Poo, see the 
dashed line in Fig. [TJ Since Pqo becomes a step function 
when s tends to infinity, the peak of k will become higher 
and narrower and finally diverge at r\ c . 

In realistic case, the system is of finite size determined 
by the size of the colliding nuclei. The finite-size scaling 
method [ill ], investigating the scaling of quantities at r\ c 
as a function of system size, is adopted to extract values 
for the critical exponents. 

The finite-size scaling of Poo suggests [HI 



P^ 



(7/ - rj c ) s 1/v 



for large s, 77 — * r\ c 



(3) 



Denoting X = (77- r] c )s l / v , $(X) is some function of X. 
v is the critical exponent of correlation length £, i.e. 



i cx 1 77 - 77 c I for r) — * r] c . 



Then, the susceptibility of P^ is 



d$ 



, v 9Poo(t7,s) l/v 

K ^ s) = —^r~ = s dx 

When ?] = rj c , ^\x=a is a constant and 
k(t] c , s) cx s 1 ^ for large s. 



(6) 



(7) 



(8) 



So the divergence behavior of k near 7? c is related with 
the critical exponent of correlation length. 

Using the critical point 7y c extracted from Pqo, we 
evaluate k(t] c ,s) as a function of s, shown in Fig. 2] 
Fitting it to Eq. © we obtain l/v = 0.739 ± 0.041, 
v = 1.353 ± 0.075. The exponent v obtained is within 
1.5% difference with the result from other calculations 
for two-dimension percolation ^fj| . 

The critical exponent of correlation length plays a spe- 
cial role in the theory of critical phenomena, because the 
scaling behavior of other quantities depends on the rel- 
ative magnitude of correlation length and system size. 
Hence, the critical exponent of correlation length ap- 
pears in the scaling relation of various quantities. By 
definition, correlation length is the distance at which cor- 
relation function reduces to 1/c [Io| . However, correla- 
tion function usually has non-monotonic behavior, e.g. in 
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FIG. 5: (Color online) D3 as a function of r\ with system size 
s = 500. Red line is a fitting by Eq.©. 



FIG. 6: (Color online) The behavior of C4 versus different 
system size. A saturation is seen. 



case of liquid, correlation function shows damping oscil- 
lation [TU , which makes the correlation length hard to be 
measured. Now we see that using the scaling behavior of 
susceptibility the critical exponent of correlation length 
can easily be obtained [l2j . 



III. APPLICATION TO RELATIVISTIC HEAVY 
ION COLLISION 

Recently, a kind of scaled third order moment of trans- 
verse momentum [T4|| 



(ptY 
(p?) 



(9) 



has been proposed as a possible signal of critical-point 
(CP) of QCD phase diagram. By imposing a tempera- 
ture gradient to the two dimensional continuum perco- 
lation of discs and assuming the transverse momentum 
for each cell takes the thermal momentum determined 
by the temperature, D^(r]) behaves like a step function, 
similar to that of P^ ual {r)) described above, cf. Fig. [5] 
and Fig. [3] This similarity makes the methods developed 
in SeclTll applicable to D3, i.e. we can locate the CP and 
determine the critical exponent of a susceptibility related 
to D 3 . 

A typical result for D3 with a half-fit to Eq. © is 
shown in Fig. [5j The inflection point C4 is regarded as the 
approximation of critical point r\ c . With the increasing 
of system size, the inflection point C4 has an asymptotic 
behavior which is shown in Fig. [SJ The saturation value 
indicates rj c = 1.124. 

Similar to Eq.Q we define the susceptibility related 
to -D3 as 



K,'(rj,s) 



drj 



(10) 



and a critical exponent v' for this thermal system as 



(11) 



Investigating the finite-size scaling of k' , 1/V = 0.642 is 
obtained, cf. Fig. [7] 
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FIG. 7: (Color online) The distribution of susceptibility k! at 
critical point r\ c as a function of system size s. The solid line 
is a fitting to Eq. (fTTjl . 



IV. CONCLUSION 

We study the critical phenomena in continuum perco- 
lation of discs. It is shown that using the inflection point 
of Poo or P£f ua as an estimation for percolation thresh- 
old is a good approximation. The critical exponent v is 
extracted from the distribution of susceptibility by using 
finite-size scaling method. 

For heavy ion collisions, different colliding nuclei have 
been tried, e.g. Pb, Au, Cu, S, C etc., which play the 
role of system size in percolation. The colliding energy 
plays the role of control variable. P^ in percolation cor- 
responds to the probability of events with QGP in heavy 
ion collision experiment, denoted as Pqgp- Studying the 
inflection point of Pqgp a s a function of y / ijv r /v, the crit- 
ical point can be extracted. A susceptibility defined as 
the derivative of Pqgp determines the critical exponent 
v which is a rare case that can experimentally measure 
the critical exponent in heavy ion collisions. 

The analogue between percolation and deconfincment 
is reasonable. The approach studied here is worth try- 
ing in the future study of the critical phenomena in real 
relativistic heavy ion collisions. 

Acknowledgement This work is supported by NSFC 
under projects No.10775056, 10835005 and 10847131. 



5 



[1] For a recent survey, see e.g., F. Karsch, J. Phys. Conf. 
Ser. 46, 122 (2006). 

[2] K. Fukushima, talk given in the seminar ol "The QCD 
Critical Point", Washington, July 28 - August 22, 2008. 

[3] R. Pisarski and F. Wilczek, Phys. Rev. D 29, 338 (1984); 
K. Rajagopal, Quark-Gluon Plasma 2, edited by R. Hwa 
(World Scientific, Singapore, 1995), p. 484; J. Berges 
and K. Rajagopa l, [hep-ph/9804233| M.A. Halasz et al, 
hep-ph/9804290; For a review, see I. Lawrie and S. 
Sarbach, in Phase Transitions and Critical Phenomena, 
edited by C. Domb and J. Lebowitz (Academic Press, 
New York, 1984), Vol. 9, p. 1. 

[4] J. Adams et al. (STAR Collaboration), Nucl. Phys. 
A757, 102 (2005); K. Adcox et al. (PHENIX Collabo- 
ration), Nucl. Phys. A757, 184 (2005). 

[5] H. Satz, Nucl. Phys. A642, cl30 (1998); H. Satz, Int. J. 
Mod. Phys. A 21, 672 (2006). 

[6] B.I. Shklovskii and A.L. Efros, Electronic Properties of 
Doped Semiconductors, Springer, New York, 1984. 

[7] A. A. Abrikosov, Spin glasses with short range interac- 



tion, Adv. Phys. 29, 869 (1980). 
[8] S. Fortunato and H. Satz, Phys. Lett. B509, 189 (2001); 

Nucl. Phys. B598 601 (2001). 
[9] K. Binder, D.W. Heermann, Monte Carlo Simulations in 
Statistical Physics, Springer- Verlag, 1988, p. 40. 

[10] K. Christensen and N. R. Moloney, Complexity and Crit- 
icality, Imperical College Press, London, 2005. 

[11] M. Barber, Finite size scaling, in C. Domb and J. 
Lebowitz, editors (1993), Phase Transitions and Critical 
Phenomena, volume 8, Academic Press, London. 

[12] D. Staufler, A. Aharony, Introduction to percolation the- 
ory, Taylor & Francis, 2003. 

[13] N.H. March and M.P. Tosi, Introduction to liquid state 
physics, World Scientific Publishing Co., Singapore, 
2002. 

[14] Du Jiaxin, Ke Hongwei, Xu Mingmei and Liu Lianshou, 
The characterizing variable for critical point in momen- 
tum space, arXiv:0810.1989 nucl-th]. 



